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Abstract 

This paper gives sharp rates of convergence for natural versions of the Metropolis algorithm 
for sampling from the uniform distribution on a convex polytope. The singular proposal dis- 
tribution, based on a walk moving locally in one of a fixed, finite set of directions, needs some 
new tools. We get useful bounds on the spectrum and eigenfunctions using Nash and Weyl- 
type inequalities. The top eigenvalues of the Markov chain are closely related to the Neuman 
eigenvalues of the polytope for a novel Laplacian. 

1 Introduction 

1.1 Overview 

The Metropolis algorithm and the Gibbs sampler (also known as Glauber dynamics) are often used 
together as one of the basic tools of scientific computation. We treat the following example: let f2 
be a polyhedral convex set in d dimensions. To sample from the uniform distribution on f2, from 
a point x in Q, pick a direction e from a fixed finite collection. Set y = x + ue where u is chosen 
uniformly in [— h, h]. If y € O, move to y. Else, stay at x. Under a mild generality condition 
on the set of directions in relation to f2, this Markov chain converges to the uniform distribution 
on Q. Our main result gives a sharp determination of the exponential rate of convergence of this 
algorithm. It is ce~ n9 ^ with g{h) asymptotic to h?v for v the first non zero eigenvalue of a novel 
Laplacian defined on f2 with Neumann condition on the boundary. 

Sampling from a convex set is a practical problem. For example, choosing a uniformly dis- 
tributed 100 x 100 doubly stochastic matrix [I] or a uniformly distributed 100 x 100 tri-diagonal 
doubly stochastic matrix [5j. It is also a basic problem of study in theoretical computer science 
[9J, [10] . Many algorithms have been proposed and studied. A readable textbook description of the 
Gibbs sampler is in Liu [8]. See [6] for a review of rigorous results for the Metropolis algorithm 
in finite spaces. The popular hit and run algorithm [H [TTJ [15] was introduced for this purpose. 
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Hit and run makes long moves and will probably be preferred in practice to the local algorithms 
studied here. 

Spectral techniques for analysis of the Metropolis algorithm on continuous spaces are developed 
in (3j El [7J . The proposal distributions there are "ball walks" choosing from the uniform distribution 
on the interior of a ball. The discrete set of directions studied here is widely used in practice and 
necessitates new ideas. Present problems can also be studied by Harris recurrence techniques 
|X2^ [T3] and by the path techniques of Yuen [16] . These give useful results but do not get the sharp 
rates on the exponents derived here. 

The remainder of this section gives a careful description of the Markov chain and the geometric 



connection between the underlying directions and the convex set Q required for ergodicity. Section 2 
gives bounds on the spectrum and eigenvectors using Nash inequalities and Weyl-type inequalities. 
Section "3| uses this spectral information to get rates of convergence. Section 4 proves that our 



operator (suitably rescaled) converges, in the strong resolvent sense, to a novel Laplace operator on 
O with Neumann boundary conditions. A similar convergence of the ball walk Metropolis operator 
to the usual Neumann Laplacian is a key ingredient of [U [7J . The final section shows how to modify 
the argument to handle a continuous choice of direction. 

1.2 Basic definitions 

Let f2 be an open convex polytope in d > 2. Thus there exists linear forms £j : M. d — > M, j = 
1, . . . , m and real numbers bj such that 

n = jx G M d , Vj = 1, . . . ,m, £j(x) > bj\ (1.1) 

Assume also that 0, is bounded and non empty. 

Consider (£ = {ei, . . . , e p } a family of vectors in M. d . For any j E {1, . . . ,p} we introduce the 
operator acting on continuous functions Mj^f(x) = mj^{x)f{x) + Kj^f(x), where 

Kj,h{f)(x) = \ f l n (x + ht ej )f(x + ht ej ) dt (1.2) 
2 J tel-l,i\ 

and rrij,h(x) = 1 - K j)h (l)(x). 

The local Metropolis operator associated to the family (£ is 

1 p 

M h (f)(x) = p J2 ( K i,hf( x ) + m jih (x)f(x)) . (1.3) 
In the sequel, denote rrih = ^ YTj=i 

m jih and K h = \ Y?j=\ K i,h- Let M h (x,dy) be the Markov 
kernel associated to this operator. This defines a bounded self-adjoint operator on L 2 (U). Moreover, 
since Mh(l) = 1, || \\l 2 ->l 2= 1- Thus the probability measure vo ^^ on ft is stationary. For 
n > 1, denote by Mj£(x, dy) the kernel of the iterated operator (Mh) n . For any x G M%(x, dy) is 
a probability measure on f2, and our main goal is to get some estimates on the rate of convergence, 
when n — > +oo, of the probability Mfi(x, dy) toward the stationary probability vo ^ n ^ ■ 

A good example to keep in mind is the case where Q = An is the set of iV x N doubly stochastic 
matrices. In other words, 



A 



N 



I (ai,j)i<i,j<N, Vi, j, ay > 0, ^a ifc = ^a fc j = 1 > . (1.4) 

I k k ) 
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The set An can be viewed as convex open polytope in A N = {{ctij)i<ij<N, Ylk a ^ = J2k a kj = !}• 
A good way to sample from An is to use the Metropolis strategy in the following manner. Starting 
from a matrix A G An choose two distinct rows Ri 17 Ri 2 and two distinct columns Cj 17 Cj 2 at 
random. Denoted = (ii, 12, jl, 32) and F = F(i) the matrix such that F{j = Si 1 j 1 —5i 1 j 2 —5i 2 j 1 +5i 2 j 2 . 
For h > given, build the family of matrices (A(t) = A + tFii))^!^^. For any t£i the matrix 
A(t) belongs to the set A N . Taking t G [—h, h] at random and keeping the move A — > A(t) only 
if it results in an element of An, we are exactly in the above situation with (£ = {F(i)}. This 
algorithm is used in [1J to study things like the distribution of typical entries or the eigenvalues of 
random doubly stochastic matrices. 

Let us go back to the general problem. From the definition of Q, a point x G M rf belongs to dfl 
iff there exists a partition I U J = { 1 , . . . , to} such that I ^ and 

Vi G I, £i(x) = h and Vj G J, ^-(x) > bj. (1.5) 

Define the following function c : R d — > N U {+00} by 

c(x) = if x G O 

= +00 ifxGM d \n (1.6) 
= card(J) if x G 90. 

To proceed, the following geometric condition is needed; it shows how the generating set (£ must 



be related to the convex set VL. Proposition 1.5 shows the condition is equivalent to having a 
spectral gap. 

Definition 1.1. The family (£ is weakly incoming to the set VL if for any point xq G d£l there exists 



e > 0, 6 G {±1} and e G (£ such that, for c defined in (1.6) 



c(x + Ote) < c(x ) ViG]0;e]. (1.7) 



The following observation is simple and fundamental. Suppose that (£ is weakly incoming, then 
span(<£) = Mr. Indeed, otherwise there is a hyperplane .ff = (IRz^)- 1 - of M. d such that span((£) C H. 
Since £1 is compact, the function x G $1 1— >■ (x, v) would have a global minimum in some xo G dVt. 
Since $7 is open, 11 C 10 + where if + = {y G (y, z/) > 0}. As (£ is weakly incoming, there 
is u G span((£) such that c(xo + u) = 0. In other words, Xo + u G D (xo + i?). This contradicts 
C x + # + . 

Example 1.1. Consider f2 t/ie convex hull of on equilateral triangle (ABC) in M? and (£ = {ei, 62} 



like on Figure 1. For a g]0, 7t/3], (£ is weakly incoming to f2 whereas for a G]7r/3, 7r[, condition 



(1.7) is satisfied in every point xq 0/ f/ie boundary excepted in point A. 



Remark 1.2. In £/ie above case of doubly stochastic matrices, the set <B = {F(i)} is weakly incom- 
ing. Indeed, if A is in the boundary of An, there exists i\ and ji such that Ai 1 j 1 = 0. Since A is 
doubly stochastic, there exists 12,32 such that Ai ± j 2 > and Ai 2 j 1 > 0. Let e = min(Ai 1 j 2 , Ai 2 jA /2, 
then for all t G]0,e], c(A + tF(i 1 ,i 2 ,ji, 32)) < c (^)- 

Denote Hk = ker(^) and let v\. be the unit vector such that £k(vk) > and H7 = {y G 
R d ,£ k (y)>b k }. 
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Figure 1: Weakly incoming condition in the case of an equilateral triangle. 

Definition 1.2. Let u G W 1 \ {0}. The vector u is incoming to H}~ if (u,Vk) > 0. Further, u is 
strictly incoming to 77^ if (u, fk) > 0; u is strictly outgoing to 77^ if (u, v^) < 0; u is parallel to 77& 
if u G -fffc- 

Lemma 1.3. Suppose that £ is weakly incoming to Q and let xq £ fl and k = c(xq). There exists 
r > and 7 C {1, . . . , m} such that |}7 = k and for B(xo, r) the open ball of radius r about x$, 

Q n B(x , r) = B(x , r) n (n ie/ 77+) . (1.8) 

Further, there exists ... , G {1, . . . ,p}, a family (0 n )n=i,...,k of numbers in {±1} and a bisec- 
tion {1, . . . , k} 3 n h-> i n G 7 suc/i i/iai /or all n 6 {1, ... , A;}, 

is strictly incoming to Hi n ; ^ 
6n£i3 n is incoming to Hi m , Vm > n. 



Proof. Proceed by induction on k = c(xo). When k = 0, there is nothing to prove. 

Suppose that the property holds true at rank k' < k — 1 and let xo G be such that c(xo) = k. 
By definition of S7, there exists r > and 7 C {1, . . . , m} with #1 = k such that 

O n B(x , r) = B{x ,r) n (n ie/ 77+) (1.10) 

Since G; is weakly incoming to S7, there exists q G {l,...,fe}, 6>i = ±1, /3i G {l,...,p} and 
ii,...,i q G 7 such that Q\ep x is strictly incoming to 77j n for n = l,...,q and is parallel 

to 77j for i G 7' := 7 \ . . . By definition of $7 there exists x' Q close to xo and r' > 

such that 7?(xQ,r') n f2 = 7?(xg,r') n (Djg/'T/^ - ). From the induction hypothesis, there exists 
Pq+i, . . . , G {1, . . . #<j+i, • • • , = ±1 and a bijection {q + 1, . . . , k} 3 n i-> i n g 7' such that 
for all n > g + 1 , 

is strictly incoming to 77 in ; 
^nC/3 n is incoming to 77j m Vm > n. 
and the proof is complete. □ 

Corollary 1.4. Suppose that (£ is weakly incoming to Q. There exist r > and e g]0, 1], such that 
for all xq G there exists q G {1, . . . ,p} and 6 q = ±1 such that 

x + t9 q e q eQ Vx G B(x ,r) nil, V* G [0,e]. (1.12) 
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Proof. The fact that r, e > can be chosen uniformly with respect to xq follows easily from 
compactness of f2. The statement is trivial when xq G f2. Suppose that xq G dQ. From Lemma 



1.3 there exists r > and / = . . . ,ik} C {1, . . . ,p} such that 

n n B(x , 2r) = B{x , 2r) D (n ie /ff+) (1.13) 



and 6*i = ±1, /3i G {!,...,£>} such that 



Sie^j is strictly incoming to H^; 
Q\e$ x is incoming to Hi q , Vg > 1. 



(1.14) 



Let x G -B(xo, r) n f2 and e G]0, r[. Then 

(fi, fiA) > 0, Vi G I => x + Wie^ G 0. (1.15) 



Thanks to (1.14), the left hand side of the above property is satisfied and the proof is complete. □ 
Proposition 1.5. The family <£ is weakly incoming to £1 iff 1 is not in the essential spectrum of 

Proof. If £ is weakly incoming to Q, 1 is not in the essential spectrum of Mh thanks to Proposition 



2.2 of this paper and Theorem 1.1 in [I]. 

Suppose now that (£ is not weakly incoming to f2. This means that there exists xo G <9f2 such 
that (1.7) does not hold. Let k = c(xo). There exists a neighborhood V of xo and I C {1, . . . , m} 
with ft/ = k such that V D = V D (D^iHi). Then, for any 8 = ±1 and any j G {1, . . . ,p}, the 
following holds true: 

If 8ej is strictly incoming to one of the (f/j)j e j, 
then Oej is strictly outgoing to one of the {Hi)i & i. 

Otherwise, there is j G {l,...,p} and 9 = ±1 such that 6ej is strictly incoming to one of the 
(Hi)i£j and incoming to the other. Then for t > small enough, c(xo + Otej) < c(xq). 
Hence, assume that there exists r > 1 such that 

• for any j G {1, . . . , r}, ej and — ej are strictly outgoing to some of the (ilj)j 6 j; 

• for any j G {r + 1, . . . ej is parallel to the (iJi)j e j. 

Recall that denotes the unit incoming orthogonal vector to Hi. Let W = span(z/j,z G I) and 

near xn use the variable x = xn + fx', x") with x' <E W and x" G W- 1 . Let y(x') = Hi , , and for 

v y v ' |<|cc'|<l 

X,h > denote / Aft (x) = (A/i)-* m W/ 2 x(^). Since any v G W x is parallel to the (fli) i6 j, there 
exists Ao, Co > such that for all h G]0, 1] and A G]0, Ao], ||/Ah||L 2 (f2) — c o- 

For any j G {r + 1, . . . is parallel to the (iTj)j 6 /. Hence, the function i i— >• f\h{x + /rfej) 

is constant and (Mj^ — l)f\h{x) = 0. 

On the other hand, for any j G {1, . . . , r} there exists G / such that ej is strictly outgoing 
to Hij and — ej is strictly outgoing to Hy . Consequently, there exists 7j, 5j > such that for t > 0, 

x G f2 and x — tej G ==> dist(x — iej, JEfy.) < dist(x, ify. ) — 7ji 

J J (1.16) 

x G f2 and x + iej G ft => dist(x + iej, JET*.) < dist(x, flj.) — <5jt. 
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Let us compute the potential mj^ on the support of f\h- For x G supp(f\h), \xj\ < Xh for all 
j = 1, . . . , r. In particular dist(a;,i?j.) < Xh and dist(a;, H^,) < Xh and thanks to (1.16), 



1 - m j h (x) 



/ ln(x + htej) + 1q(x — htej) dt 
Jo 



< I dt+ [ dt<x(— + -*r- 

lo<t<dist(x,Hi.)/(8jh) JO<t<dist(x,H i ,)/('y j h) \7j °j 



J 3 



Finally, 

((1 - M h )h h , fx h ) LHn) = li>2((l- M jth )f Xh , f. 



1 2=1 



< l jZl { l - m i,h{x))\fx h (x)\ 2 dx<CX\\f Xh \\ 2 L2 



(ay 



(1.17) 



(1.18) 



Here we used the fact that for any non-negative fonction /, one has {Kj^f, f) > 0. Finally, we 
conclude by taking A = 2~ n — > as n — > oo. Indeed, the functions ji-^h are mutually orthogonal. 
Their norm is bounded uniformly from below and they satisfy < ((1 — M/ l )/ 2 -n ft , / 2 -™h) < 
C2~ n . □ 



2 Spectral Analysis of the Metropolis Operator 

This section is devoted to the analysis of the spectral theory of the Metropolis operator. For this 
purpose, we introduce a Laplace operator associated to the family <£ to be used as a model. For 
any e G M. d \ {0} and any smooth function u, define d e u(x) = ^(u(x + te))\ t= Q. Then, consider the 
operator Ag, defined by 

1 P 

3 (2.1) 
£>(A g ) = {ue H^fi), A^u G L 2 ,d n ^u m = 0} 

with d n ^u(x) = i X]j=i( n ( x )' e j)de u{x), n(x) denoting the outgoing normal vector to the bound- 
ary at point x. If the domain Q has smooth boundary, the normal derivative is well defined. In the 
case where it is Lipschitz, it can be defined by duality in the following way. 

Define first the gradient and divergence associated to the family (£, by divgit = ^X^j=i^ej M i 
for any u = (u%, . . . , u p ) and V^u = (d ei u, . . . , d ep u). Then, define a trace operator 7u by 

7l£ : {u G {L 2 {n)) p , div E (u) G L 2 (ft)} -> H^ 1 ^ 2 (dQ) (2.2) 
and and for v G H 1 ^), 

div(j(n)(x)f (x) dx = — / (u(x),\7 t ev(x)) <CP dx+ j ^^{u)v\Q^da{x). (2-3) 
P Jn Jdn 

In particular, for u G H 1 ^) satisfying A^u = gdiv^Vgu G L 2 (S7) define d n ^u\QQ = 7u(V(cu) G 
H^ 1 ^ 2 (dQ) and the set .D(— Ag) is well defined. The Dirichlet form associated with — Ag is 

1 p r 

W = r V / \d ej u{x)\ 2 dx. (2.4) 
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Let <£q be the canonical basis in IR d . Then, A<j = where A is the usual Laplace operator 
and £{£ {f) = id In \^ f\ 2 dx is the usual Dirichlet form. Since <£ = {ei, . . . ,e p } spans a simple 
calculation shows that there exists a constants C > such that 

C-^fiKSeWKCS^f). (2.5) 

Then, it is standard to show that — A<j is the self-adjoint realization of the Dirichlet form £<t. 
A standard argument using Sobolev embedding shows that — A<e has compact resolvant. Denote 
its spectrum by vq = < v\ < i>i < . . . and by mj the associated multiplicities. Observe that 



mo = 1. Section 4 shows that h (1 — M^) converges to — Ag in the strong resolvent sense so that 
eigenvalues and eigenvectors converge; see [13] . 
The main theorem of this section follows. 

Theorem 2.1. Suppose that (B is weakly incoming to £1, then the following hold true. 

i) There exists ho > 0, 5q G]0, |[ and a positive constant C such that for any h e]0, ho], the 
spectrum of is a subset of [— 1 + So, 1], 1 is a simple eigenvalue and Spec(M^) n [1 — 5q, 1] 
is discrete. 

ii) For any h £}0,ho] and < A < 5oh~ 2 , the number of eigenvalues of in [1 — h 2 X, 1] (with 
multiplicity) is bounded by C(l + \) d / 2 . 

Hi) For any R > and e > such that Vj+i — v-j > 2e for Vj+2 < R, there exists h\ > such that 
one has for all h s]0, h\], 

Spec ( - ~ h f h ^j n]0, R] C Uj>i[vj - e, Vj + e], (2.6) 

and the number of eigenvalues of X ~^ h in the interval [vj — e, Uj + e] is equal to rrij. 

A consequence of this theorem is that Mh has a spectral gap g(h) = 1 — sup(Spec(M^)\{l}) > 
and that lim^_ s>0 + a 2 9(h) = v \- This will be used in the proof of total variation estimates. 



The strategy used to prove the first part of Theorem 2.1 is very close to the one given in [3]. 
First, show that some iterate of the Markov kernel "controls" the random walk on a ball. Next, 
this ball walk on the polytope is compared to the same walk on a large torus containing fi. Finally 
the information on the torus is transferred back to the original problem. 



The proof of the last part of Theorem 2.1 is slightly different from the proof in [3]. Indeed, the 
starting point of the analysis in jJJ is that for regular function ip with normal derivatives vanishing 
on the boundary, h~ 2 (l — Th)ip is close to —Aip up to the boundary, where Th is the Metropolis 
operator associated to the kernel vol(.B(0, l))~ 1 h~ d l\ x _ y \ <h . Here, this property fails to be true. 
Suppose for instance that SJcl 2 and that its boundary is given near (0, 0) by x\ > 0. Suppose 
that e\ = (a, b) and ei = (6, —a) for some a, b > 0. Then 

h- 2 (l - M 1 . h )f(x) = / (f(x) - f(x + tei)) dt 

Zn J \t\<h,x+te!eQ 

1 d e J(x) [ tdt + 0(l) (2.7) 



2/i 3 



|t|</i,a!+£eien 

2' 



^deJ(x)l ]0 , ah] ( Xl )(h 2 -^j +0(1). 
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A similar expression holds for M2 h an d summing these equalities gives 



h-\l - M h )f(x) 



1 



deJ(0,X2)l]0 y ah](xi) [h? - - 

+ d e2 f(0,x 2 )l] ,bh](xi) I h 2 



b 2 



+ 0(1) (2i 



If a = b, d ei f + d e2 f is proportional to the normal derivative of / and hence, the above quantity 
is bounded. 

Suppose now that a < b. Then the above quantity is bounded on x\ £ [ah, bh] provided 
d e2 f(0,X2) = 0. Then the same argument on [0, ah] shows that d ei f(0, xz) = also. 

In order to avoid these difficulties, we work directly on the quadratic form and show that the 
Dirichlet form associated to the Metropolis operator converges to the Dirichlet form of the Laplace 
operator with Neuman boundary conditions. The end of this section is devoted to the proof of 
Theorem 12. 1[ 



(2.9) 



Proposition 2.2. There exists N £ N and constants c\,C2 > such that for all h e]0, 1] 

M h ( x , dy) = fih(x, dy) + cih~ d l\ x _ y \ <C2h dy 
where for all x £ f2, fj,h(x,dy) is a positive Borel measure. 

Proof. The proof follows the lines of |4]. Denote = - YTj=i Kj,h- Since for any hi > hi > and 
any non-negative function /, h-zK-hiJ > hiK^f, it is sufficient to prove the following: there exists 
ho > 0, ci, C2 > and N £ N* such that for all h e]0, ho], one has, for all non-negative continuous 
functions /, 



K?(f)(x) > Cl h~ d [ 

Jy 



y£Q,\x—y\<C2h 



f(y) dy. 



(2.10) 



First note that it is sufficient to prove the weaker version: for all x° G Q, there exist N(x°),a = 
a(x ) > 0, c\ = ci(xo) > 0, C2 = C2(xo) > 0, ho = ho(xo) > such that for all h g]0, ho], all x £ 17 
and all non-negative functions / 



x 



x u \ < 2a 



< ( * o) (/)(x) > Cl h 



-d 



j/£f2,|x— y\<C2h 



f(y)dy. 



(2.11) 



Let us verify that (2.11) implies (2.10). Decreasing a(xo) if necessary, it may be assumed that 
2a(xo) < r(xo), where r(xo) is given by Lemma 



1.3 



Since Q is compact, there exists a finite set 
F such that Q, C U Xq ^f{\x — xo\ < a(xo)}. Let ./V = sup{A r (xo), xo £ F}, c\ = mm XoG p Cj(xo) 
and h' = min XQ& F ho(xo). One has to check that for any xo £ F and any x with \x — xo\ < a(xo), 
the right inequality in (2.11) holds true with N = N(xq) + n in place of N{xq) for some constants 



c\,C2,ho- Moreover, one may assume that /iomax|ej| < min a . ( 
Let e > 0, q £ {1, . . . ,p} and 8 q = ±1 be given by Corollary 
one has 



(z F u(x )/N. 

Then for |x—xq| < (2— jj)a(xo), 



1.4 



K 



N{x )+l 



1 



f(x) > p K hA K h 



N(x ) 



PJo 



K 



N(x ) 



P 

> c' h- d 



f(x + ht9 q e/3 q ) dt 

f(y) dydt 



y£fl, \y-x-ht8 q ej3 q \<c 2 h 



(2.12) 



yen, \y-x\<c 2 h/2 



f{y) dy 
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since for any t £ [0,min(e, 2ma ^ 2 j e | )], {\y — x\ < C2/1/2} C {\y — x — ht9 q ep q \ < czh}. Iterating this 
computation n < N times gives (2.10). 

It remains to prove (2.11). If xo £ Q, the proof is obvious. Indeed, since (£ spans IR^, it is easy 
to see that for any 5 > 0, there exists C3, C4 > such that for any non- negative function /, 



dist(y, an) > Sh => Ki(f)(y) > c 3 h- d [ f(z) dz Vy G fi. 



(2.13) 



Suppose that xq G an and denote = c(xq). Let (ij)i<j<k, {fij, 



i<j<ki \ a j)i<j< 



k be as in Lemma 



1.3. Let 1 = 71 > 72 > • • • > 7& > and 61, ... ,6k > be such that for all j, jj — 6j > 77+1. Let 



Gj = [7j — 6j,jj] and G = H k =l Gj. In the following computation, c denotes a positive constant 



3= 

independant of / and h that may change from line to line. Since / is non-negative 

k 

1 f\ ( _\ ^ „_— k 



K(f)(x) > p- k K^ h . . . K^ h f(x) > c 



fix + h^djtjefy) dt 

t&A h {x) j=1 



(2.14) 



where A h {x) = {t = (ti, . . . , i^) € G, Vi = 1, . . . , k, x + h Y,j=i Qjtj e Pj G 

Since B\ep x is strictly incoming to there exists some constant c§,cq > such that for any 

t e I, 



dist \x + h Qjtje^ , > c 5 /?ii - c 6 /i(t 2 H h t k ) 

3=1 J 

> c 5 /i(7i - 61) - c 6 /i(72 H h 7fc) 

> c 5 /i( 7 i - 6x)/2 



(2.15) 



by taking 72, • • • ,7fc small with respect to 71. Similarly, by taking jj very small with respect to 
Jj+i for j = 2, . . . , k, there is cj > such that for any j = 1, . . . , k, 



V(ii, ...,%) € Gi x • • • x Gj, dist I x + OiUe^,^ \n j > c 7 /i. 



(2.16) 



Hence, 



and for any TV > 



K k h f( x ) >cJ^jL + hJ2 Qjtm j dt 

K k+N f{x)>cf^K^(f) (x + hJ^O&e^ dt. 



K d h +k f(x) > ch- d [ 
Jte 



teG J \x+h Yj] =1 tjep. -y\<cg,h 



f{y) dydt. 



(2.17) 



(2.18) 



Combining (2.13), (2.16) and (2.18), there is cs > small enough such that any y £ M such that 
\x + hY^j=i tj£/3j — y\ < cg/i belongs to n and hence 



(2.19) 
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Since, K%f(y) > p k K hA . . . K hA f(y), then 



K d h +2k f(x) > ch- d [ 
J(t 



where 



B h (x) 



f(y - s i e Pj ) dtdsd v 

(t,s,y)&B h (x) j=l 



(t,s,y) 6 G X G X R d , \x + h V] Ojtje-Pj — y\ < csh and 



Using the new variable z = y — h Y^,j=i @j s j e Pj > 



/(z) dtdsdz 



(t,s,z)£D h (x) 



with 



(t,s,z)eGxGxSl, |x + /t y^(tj — Sj)9je^ j — z\ < c$h and 



j'=i 



Since in the above integral, — Sj-| < <5j, taking the Sj's small enough gives 

r 

Dh{x) D < (i,s,z) £ G x Gx < c%h/2, VI = 1, . . . ,k, z + h ^ OjSjZpj 6 ^ 

I i=i 



Now using (2.16), it follows that 

-Dh(x) D {(i, s,z) £ Gx Gx < c 8 /i/ 2 } . 



Combined with (2.22), this yields the announced result. 



(2.20) 



VZ = 1, . . . , k, y - h^2 e i s 3 e Pj e n \ • ( 2 - 21 ) 

3=1 



(2.22) 



VZ = 1, . . . , k, z + e 3 s i e Pj G n \ ■ ( 2 - 23 ) 



(2.24) 



(2.25) 
□ 



Following the strategy of [I], introduce the Dirichlet form associated to the iterated kernel M^: 



1-Mhu,u 



(2.26) 



Also, put Q in a large box B =] — A/2,A/2[ d and define an extension map E 1 : L 2 (S1) — > L 2 {B) 
which is continuous from H l {yi) into H l (B) and vanishes far from $7. This is possible since dCl has 
Lipschitz regularity. Finally, introduce the Dirichlet form on B: 



£ h (u) = h~ 



BxB,\x-y\<h 



\u(x) — u(y)\ 2 dxdy. 



(2.27) 



Then Proposition 2.2 easily yields the following (see [I] for details). 
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Lemma 2.3. There exists Co, ho > such that for any h g]0, ho] and any u G L 2 (VL), 

4 (E(u)) < C (£ htN {u) + h 2 \\u\\l Hn ^j . 

Moreover, any function u G L 2 (Cl) such that 

W u \\h(n) + h ~ 2 ((! - M ^) n ' M ) L2(n) < 1 
admits a decomposition u = ul + G ii -1 ^), < C\, and \\uh\\l 2 ^ Ci/i. 



We are now in position to prove the first part of Theorem 2.1 First, assume that M^u = u. 
Then, it follows from Proposition 2.2, that 

-d f t„.r„\ „./.^2 j„j„ / / f„.f„\ „.f„.\\2 tutN 



cih~ d / (u(x) - u{y)y dxdy < / (u(x) - u(y)) z M£ (x, dy) dx. (2.29) 

J QxQ,\x—y\<C2h JQxQ 

On the other hand, the right hand side in the above inequality is equal to Shj^[u) which is actually 
equal to zero. Hence, u is constant and 1 is a simple eigenvalue. 

Using the Markov property of , positivity of \Xh and the fact that dSl has Lipschitz regularity, 
easily yields 

W^hh^^L^ = fJ>h(fl) < 1 - ci/i" d min / l\ x - y \ <C2h dy < 1 - S' (2.30) 

xsn Jn 

for some S' > independent of h. Working as in the proof of Theorem 1 in [3] shows that there 
exists 5q G]0, |[ such that for any u G L 2 (Q) and any n> N, 

{MZu,u) L2{n) > (-l + 5o)\\u\\ 2 L2{n) . (2.31) 

Hence, the same holds true for n = 1 with a possibly different <5o- 

To show that there is 5q > sufficiently small so that the spectrum of Mh is discrete in [1 — 5o, 1] 



it suffices to work as in the proof of Theorem 4.6 in [3], using again Proposition 2.2 



Similarly, the Weyl bound on the number of eigenvalues follows from Lemma 2.3 as in Lemma 
4.8 in [3]. This proves Part i. 

To prove the last part of the theorem, work on the Dirichlet form is needed. In the following, 
denote Eh = £h,l- Introduce the bilinear form associated with Eh- 

B h (u,v) = ((1 - M h )u,v) L2(n) , Vu, v G L 2 (ft). (2.32) 
A standard computation shows that Bh(u, v) = ^ Sj=i Bj,h(u, v) with 



Bj^hiu, v) = —r f (u(x) — u(x + tej)) (v(x) — v(x + tej)) dxdt (2.33) 

4" ■'xen,x+te J en ,\t\<h 



Lemma 2.4. Let 9 G C°°(r2) be fixed and let (iph,fh) G H 1 ^) x L 2 {£1) be such that \\rh\\L 2 {Q) 
0{h) and tph converges weakly in H x (p) to some (p. Then 



and 



lim h- 2 B h (r h ,0) = (2.34) 
lim h- 2 B h (ip h ,0) = ^ / <V^(x),V^(x)> CP dx. (2.35) 
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Proof. To prove (2.34), observe that since 9 is smooth, 



h- 1 r 

(1 - M jA )9(x) = — / (9(x) - 9(x + tej)) dt 

1 J\t\<h,x+t e:j en 

deO{x) f , a , (2 - 36) 

/ tdt + 0(h 2 ). 



p 4 (x) = MW/ t « 



2*1 7[t|<h,a:+te^en 

Denoting 

/l I. f tA = — 

2^ 7[t[<ft,a!+te3-en 

observe that supp( / o/ l ) C {x G $7, d(a;,<9Sl) < /i} and Hp^IIl 00 = 0(h). Hence ||p/i||l2 = 0(h 3 / 2 ) and 
since ||r/ l ||£ ) 2 = 0(h), it follows that 

h- 2 B jA (r h ,6) = h- 2 {r h ,(l-M j>h )9) L 2 = (h^r^h' 1 p h ) L * + 0(h) = 0(h 1 ' 2 ) (2.37) 

which goes to zero as h goes to zero. 
To prove (2.35) observe that 

B(x + tej) -9(x) = tip(t,x) 

f 1 (2-38) 

tp h (x + tej) - ifh(x) = t / d ej <p h (x + tzej) dz 

Jo 

with ip(t,x) smooth and ip(0,x) = d ej 9(x). Hence 

h~ 2 Bj. h (ip h ,9) = — -g / t 2 d e .(p h (x + tzej)ip(t,x)dtdzdx 



1 



xeU,x+tejEU,\t\<h,ze[0,l] 

2 ; 



■u d ej (fh(x + huzej)tp(hu,x) dudzdx (2.39) 

4 Ja;efi,a;+/iuejen,|u|<l,ze[0,l] 



u 2 d ej fh(x)^(hu, x — huzej) dudzdx. 

x— huzej&Q,x+hu(l— z)ej&fl,\u\ <l,z(L [0,1] 



Taylor expansion of tp shows that ip(hu,x — huzej) = d e .9(x) + 0(h). Hence, for any 5 > and 
any/iG]0,l], 



h 2 Bj^ h (ip h ,9) = - I u 2 d ej ip h (x)d e .9(x) dudzdx + 0(h) 
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x— huzej £Q,x+hu(l— z)ej£fl,\u\ <l,.z6[0,l] 



= I s (h) + J 5 (h) + 0(h) (2.40) 

with /<5(/i) equal to the above integral over d(x,dft) > 5 and J$(h) the integral over d(x,d£l) < 5. 
Then, by Cauchy-Schwartz, |J,s(/i)| < C(9)5 1 ^ 2 \\ifh\\H 1 - On the other hand, for any h €.]0,S[, 

1 



6 



Is(h) = - / d e .iph(x)d e .9(x) dx 



x£Q,d(x,dQ)>S ^ 



\ I d e .<p h (x)d e .9(x) dx + O f<y 1/2 |bh|| H i 



6 



xefi 
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Given e > 0, it is easy to find 5 > small enough such that for any h G]0, S[, \Js(h)\ < e and 
\h(h) — g J^gQ d ej iph(x)d e:j 0(x)dx\ < e. Now make h — > + , 5 being fixed, and use the fact that iph 
converges weakly in H 1 to get 



lim h 2 B j h (ip h ,9) = - [ d e .<p(x)d e .9(x) dx 



and the proof is complete. 



(2.42) 
□ 



To complete the proof of Theorem 2.1, denote |A 



h- 2 (l - M h ). Let R > be fixed and 
observe that if Vh G [0, R] and fh G satisfy lA/J/h = Vhfh and ||/h||ia = 1, then, thanks to 

<Ph + r h with ||r/i|| L 2(Q) = O(Zi) and cph bounded in 



Lemma 2.3 
H . Hence 



fh can be decomposed as fh 
(extracting a subsequence if necessary) it may be assumed that iph weakly converges in 

that for any 



2.4 



H 1 to a limit ip and that Vh converges to a limit v. It now follows from Lemma 

^-JjVef(x),V & 8(x)) CP dx = v(ip,6) L 2. (2.43) 

Since 9 is arbitrary, it follows that (— Ag — z^)<^ = and d n ^(f\gQ = 0. In fact, this also proves that 
for any e > small, there exists h e > such that for h G]0, h e ], one has 

Spec(| A h |) n [0, R] C Uj[uj - e, Vj + e] (2.44) 

and 

ttSpec(|A fc |) n [uj - e, Vj + e] < m,- (2.45) 

In fact, there is equality in ( 2.45 ) . The following proof is a simplification of the one in [4] . Proceed by 
induction on j: let e > 0, small, be given such that for < Vj < M+l, the intervals Ij = \vj—€, Vj+e] 
are disjoint. Let (Hj)j>o be the increasing sequence of eigenvalues of |A/j|, on = J2f=i m j an d 
(efc)fe>o an othonormal basis of eigenfunctions of — Ag such that for all k G {1 + <7tv, . . . , on+i}, one 
has (— Ag — z^v+i)e/c = 0. As is a simple eigenvalue of both — A<j and |A/i|, clearly uq = fio = 
and m = 1 = )JSpec(| A^l) n [v$ - e,u + e\. 

Suppose that for all n < N, m n = fjSpec(| A/J) n \v n — e,v n + e]. Then by (2.44), for /i < /i £ , 

»1 +aN > - e. (2.46) 
By the min-max principle, if G is a finite dimensional subspace of H 1 with dim(G) = 1 + trjv+i, 



0"AT+1 



< 



sup {\A h \ip,i/)) L 2( Q )- 



ip&G,\\ip\\=l 

Let G be the vector space spanned by the e k , < k < o"at+i- Then, dim(G) 
follows from Lemma 2.4, for any k, k' < 1 + cjv+i, 

1 



lim /i B h (e k ,e k >) 



Hence 



lim /i 

7i-»0 



-2 



1 

ftp 



(V c e fc (x), V £ e fc '(x)) CP dx. 



| VgV'(^| 2 dx < vn+i 



(2.47) 
1 + ctjv+i and it 

(2.48) 
(2.49) 



for any ip G G with ||-t/)||^2 = 1. Since G has finite dimension, a standard compactness argument 
shows that there exists h t > such that for any h G]0, /i e ] and any ip £ G with HV'Hl 2 < L 

/i-X^VO <^JV+i + e. (2.50) 

Therefore /^ CTJV+1 < vn+i + £• Combining this with (2.46) and (2.45) gives rrijv+i = J}Spec(| Ah\) H 
[z^jv+i — e, i/jv+i + e]. The proof of Theorem 2.1 is complete. 
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3 Total Variation Estimates 



This section gives estimates on the convergence speed of the iterated kernel M£(x,dy) towards 
its stationary measure y^bx ■ Recall that the total variation — ^Htv between two probability 
measures [i and v on f2 is defined by 



Ha* - v\\tv = sup \n(A) - v{A)\ 
where the sup is taken over all measurable sets A. Equivalently, 

K/)-K/)l- 



1 



/i - V\\TV 



sup 

/e£°MI/IU 



(3.1) 



(3.2) 



Theorem 3.1. Assume that <£ is weakly incoming. Then there exists C > and ho > such that 
for all h s]0, ho] and all n £ N, the following estimate holds true, with g(h) the spectral gap studied 
in \Section 2c 



sup 



M%{x,dy) 



dy 

voi(n) 



< Ce" n9i - h) . 



(3.3) 



TV 



Proof. The proof is very close to the proof of Theorem 4.6 in [3] and is just sketched for the reader's 
convenience. Observe first that n > h~ 2 can be assumed, since otherwise the estimate is trivial 
thanks to the lower bound on the spectral gap. 

Let IIo be the othogonal projector in L 2 (S7) on the constant functions. Observe that 



2 sup 



M%(x,dy) 



dy 



vol(O) 



TV 



\ M h - H-o\\l^^l° 



(3.4) 



Using the spectral decomposition of M^, let < Ai^ < • • • < Xj^ < • • • < h~ 2 5o be such that 
the eigenvalues of in the interval [1 — 5q, 1] are the 1 — h 2 \j^ with associated orthonormalized 
eigenfunctions Mh(ej,h) = (1 — h Xj^hjejft- 

Then write — IIo = i + 2 + Mh 3 , so that the operators 1 , 2 have kernels 



M h ,i{x,y)= ^2 (1 - h 2 X jth )e jjh (x)e jih (y) 
Ai,h<A iA </i-« 

M hj2 (x,y) = ^2 (1 - h 2 \j >h )e jih (x)e jth (y) 



(3.5) 
(3.6) 



where a G]0, 2] is a small constant that will be chosen later. Then 

3 

"\\M1 , 



2 sup 



M%(x,dy) 



dy 



TV j=1 



(3.7) 



vol(fi) 

and terms on the right hand side must be estimated. 

From (2.30), it is easy to prove that any eigenfunction Mh(u) = \u with A G]l — do, 1] satisfies 

\\u\\ L °° < Ch~ d/2 \\u\\ L 2. (3.8) 

As in (3], using in particular the bound on the number of eigenvalues, we show that for n G N, 

I|M£ 2 ||l^l~ + ||M£ 3 || L oo_ LOO < C((l - h 2 ~ a ) n + (1 - 5 ) n )h- 3d / 2 (3.9) 
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For n > h 2 , this implies that 

\\^h,2\\L°°^L°° + \\M^ 3 \\l°°-^l°° < C a e~ nh . (3.10) 

It remains to estimate Let E a denote the space spanned by the eigenvectors ej t h such that 

h — h~ a . Then, thanks to Part ii of Theorem 2.1, dim(E a ) < h~ da l 2 . As in [3], Lemma 



shows that there exists a > and p > 2 such that for any u G E a , 

\\u\\ 2 LP < Ch~ 2 (£ h>N (u) + h 2 \\u\\ 2 L2 ) . 



2.3 



This gives the following Nash estimate, with j> = 2 

i 



> 0: 



,2+" 



\u\\ L2 D < Ch- 2 (£ hjN (u) 



t 2 ii \\2 

h m 



Vu £ E a . 



(3.11) 



(3.12) 



li 2 ) H^lLi 

This inequality allows an estimate of Mft i from L 1 into L 2 and this leads to HM^H^oo^^oo < 
Ce -kN 9 (h) for k > ^-2 Ag Mh ig bounded by 1 on L°° it follows that kN can be replaced by 
n > /i~ 2 in this estimate, and the proof of Theorem 



3.1 



is complete. 



□ 



4 Convergence of the Resolvants 



Let us denote |A^| = h (1 — M^). Recall Ag from (2.1). This section proves strong resolvent 



convergence of |A/j| to Ag. For background and consequences, see |13j . 
Theorem 4.1. Let z G C \ [0, +oo[ and 5 G L 2 (ft). T/ien 



lim HflAfcl -z)" 1 5 -(-A e -z)~ 1 5 || r9 _ N =0. 



/i->0+ 



(4.1) 



Proo/. Let z G C \ [0, +00 [ and 5 G L 2 (Q) be fixed. For any h > let G L 2 (f2) be the solution 
of (|A/j| - z)/ h = 5. Hence 



-z{f h ,fh)L2 + 



1-M fe 
/i 2 



L-' 



h/L 2 - 



(4.2) 



Since z ^ [0, oof and |A/j| is a positive operator, it follows that H/hHz^ < dist(z, [0, oo[) _1 ||<7|| L 2 is 
bounded uniformly with respect to h. It follows from the above equation that there exists Cq > 
such that 

\\f h \\ 2 L 2+h- 2 £ h (f h )<C \\g\\l2. (4.3) 



It now follows from Lemma 2.4 that there exists C > depending on z and ||g||£2 such that for 
any h g]0, 1], we can write fh = <Ph + fh with H^Hjji < C and H^Hl 2 < Ch. Let / G H 1 ^!) 
and (hk)ken be a sequence of positive numbers such that ((fh k )k converges weakly to / in H . Let 
6 G C°° {U) be fixed. Then 

-z{f hh ,0) L2 +h- 2 B hh (f hk ,e) = (g,9) L2 



(4.4) 



and taking the limit k — > 00 it follows from Lemma |2.4| that 

1 



g(x)6(x) dx. 



(4.5) 



Since # is arbitrary, this implies (— Ag — z)/ 
subsequence (/ifc), this shows that \\fh — J\\l 2 ~ 



- g and d n ^f\Q^i = 0. Since, this is true for any 
when h — > 0, which is exactly (4.1). □ 
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5 Some Generalizations 



Here we present a possible generalization of the previous results. It is still assumed that is a 
convex polytope in M. d . Suppose that E C M. d is endowed with a Borel probability measure [i. For 
any e £ F, define 

K e ,hf{x) = \ I f(x + hte)dt (5.1) 

1 Jte[-l,l],x+hteeQ 

and 

K h f(x) = [ K eth f(x)dn(e). (5.2) 

JeeE 

The associated Metropolis operator is defined by M^/(x) = rrif l (x)f(x) + K/ l f(x) with irih(x) = 
1 /uil). 

Definition 5.1. Say that (E,(i) is weakly incoming to Q if for any xq £ <9f2 there exists e > 0, 
9 £ {±1} and a measurable subset F C E such that /i(F) > and 

c{x + 6te) < c(x ) Vt e]0, e], Ve £ F (5.3) 



Lemma 5.1. There exists some measurable subsets F\, .. .,F^ C E such that /^(Fj) > for all j 
and any (/i, . . . , fy) £ ILy =1 Fj spans M. d . Moreover the sets Fj can be chosen with arbitrary small 
diameters. 



Proof. From the same argument as in remark following Definition 1.1 we can easily see that \x can 
not be supported in an hyperplane of M. d . Let us prove by induction that for k = l,...,d, there exists 
Fl, . . . Ffc C E such that fi(Fj) > for all j and for any (/i, . . . , fk) £ Hj =1 Fj, rank(f±, . . . , fk) = k. 

If k = 1, it suffices to take Fl C F \ {0} with /x(Fi) > 0, which is possible thanks to the fact 
that F is weakly incoming to Q. 

Assume that the property holds true at rank k — 1 < d. There exists Fl, . . . F^i C E such that 
n(Fj) > for all j and any (fx, ... , fk-i) £ II^ZjFj, F~ = span(fi, . . . , /fc_i) has dimension k — 1. 
Since supp(^) is not contained in H, there exists F& C F\H with fi(Fk) > 0. Then Fi, . . . , F& 
satisfy the property at rank /c. 

The fact that we can take diam(Fj) arbitrary small can be shown as follows. Let e > and 
assume by contradiction that there exists jo such that for any / £ Fj , fi(B(f, e) n Fj ) = 0. Then 
any compact subset of F J0 would have measure zero, which is impossible since /Lt(F JO ) > 0. □ 

Introduce the following differential operators associated to the set F: 

V E ■ H\Q) -»■ L°° (E, L 2 (ft)) (5.4) 

defined by V^u(e,x) = (Vu(x),e) C d for any (e,x) £ F x $7; 

divij : L°° (F,iF(Q)) L 2 (ft) (5.5) 

defined by div^/(a;) = f E (V x f(e, x), e) C dd(i(e) for any x £ f2; and 

A £ : F/ 2 (0) ^L 2 (0) (5.6) 

given by A^ = gdiv^Vs. 
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Define also the following trace operator: 

1e : {/ G L°° {E,H\n)) , d\y E f G L 2 (Q)} -> fT5(Sfi) (5.7) 

by 

/ lEf( x ) v \an(x) do(x) = \ div E f(x)v(x) dx + / / f(e, x)V E v{e, x) dxd^i(e) (5.8) 
JdQ Jn JeJo. 

for any v G JT^fi). Observe that if / G L°°(E, C 1 ^)), then 

lEfi x )= I {e,n(x)) cd f(x,e)dfj,(e) (5.9) 



where n(x) denotes the unit outgoing normal vector to the boundary d£l at point x. 

For u G H 1 ^) such that A^ii G L 2 (J7), the function / = Vgw satisfies div^/ G L 2 (S7). Hence, 
the operator 

7^:{u£ fl 1 ^), A E u G L 2 (fi)} -> £T5(8n) (5.10) 

defined by 7^u(x) = 'y E V E u(x) is continuous. 

Finally, introduce the following quadratic form on i? 1 ($7) : 

£ E (u) = - [ [ \V E u(e,x)\ 2 dxdn(e) VuGH 1 ^). (5.11) 
6 JeJq 



From Lemma 5.1 it follows that, since E is weakly incoming to Q., there exists some subsets 
Fi, . . . , Fd with arbitrary small diameters and fJ>(Fj) > such that any (/i, . . . , fy) G F\ x • • • x 
spans M. d . Taking the diameter of the Fj sufficiently small, it is easy to show that there exists C > 
such that for any u G i? 1 (f2), 

^l|Vu||£ 2(n) < 8 E {u) < C\\V U \\ 2 L2{ny (5.12) 

Then, the operator —A E = — gdiv^V^; with domain D(—A E ) = {u G H 1 ^), A E u G 
L 2 (Q,), j E u = 0} is the self-adjoint realization of the Dirichlet form £ E . Moreover, it follows 
from (5.12) that — A^ has compact resolvant. Denote its spectrum by uq = < v\ < vi < . . . and 
by rrij the multiplicity associated to Vj. Observe that m,Q = 1. 

Theorem 5.2. Suppose that (E,fx) is weakly incoming to fl, then the following hold true. 

i) There exists ho > 0, So G]0, ^[ and a positive constant C such that for any h g]0, ho], the 
spectrum of Mh is a subset of [—1 + So, 1], 1 is a simple eigenvalue and Spec(M^) D [1 — So, 1] 
is discrete. 

ii) For any h G]0, ho] and < A < 5oh~ 2 , the number of eigenvalues of Mh in [1 — h 2 X, 1] (with 
multiplicity) is bounded by C(l + \) d / 2 . 

Hi) For any R > and e > such that Vj+i — Vj > 2e for Vj+2 < R, there exists h\ > such that 
one has for all h g]0, hi], 

Spec ^ n]0, R] C Uj>i[vj - e, u, + e] (5.13) 

and the number of eigenvalues of 1 ^ h A 2 Ih in the interval [vj — e, Uj + e] is equal to rrij. 
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Here are two examples of (E, fj,) which are weakly incoming to 0. The first is the case where 
E = {ei, . . . , e p } is discrete and p is simply the measure - Y3i=\ ^e=e r Then it suffices to assume 



that E is weakly incoming to f2 in the sense of Definition Moreover, in that case the conclusion 



of Theorem 5.2 are exactly those of Theorem |2.1| 

A second example is the following. Let E be equal to the sphere S"^ -1 and p = dad be the surface 
measure. Assume that p : S^ 1 — > M + is a continuous function such that f S d-i p{u)dad{u) = 1 and 
let p = p(uj)dad(uj). Then (E,p) will be weakly incoming to f2 iff there exists a family of vectors 
ei, . . . ,e p E supp(/>) such that (ei, . . . , e p ) is weakly incoming in the sense of Definition 1.1 For 
instance, if p is strictly positive on S d ~ l then these assumptions are automatically satisfied. 



The proof of Theorem 5.2 is very close to that of Theorem |2.1| and only the main steps are 



given. The following proposition is a version of Lemma 1.3 adapted to the present setting. 



Proposition 5.3. Assume that (E,p) is weakly incoming to $7, let xq £ fi and denote k = c(xo). 
There exists e > and some subsets Fi, . . . , C E such that p{Fi) > for all i = 1, . . . ,k and 

• there exists ro > and I C {1, . . . , m} with jj I = k such that 



nnB(x ,r ) 



n B(x ,r ); 



(5.14) 



• there exists 9%, ... ,9k S {±1} and a bisection {1, . . . , k} 3 n i— > i n € I such that for any 
n = 1, . . . , k and any f n E F n , 

9 n fn is strictly incoming to Hi n (5.15) 

and 

Onfn is incoming to Hi m Vm > n. (5.16) 
Moreover the sets F\, . . . ,Fk can be chosen with arbitrary small diameter. 



Proof. First, it is clear that (5.14) holds true. We prove (5.15) and (5.16) by induction on k = c{xq). 
For k = there is nothing to prove. 

Assume now that the property holds true for all x' such that c(x' ) < k — 1 and suppose that 
c(xo) = k. Since (E,p) is weakly incoming to f2, there exists F C E, 9\ G {±1} and e > such 
that c(xo + t9\f) < c(xo) for all t e]0, e]. Assume without loss of generality that 9\ = 1. Since 
p{F) > 0, there exists /° € F such that for all p > 0, p(B(f°,p) n F) > and 

c(x + t/) < c(x ) Vf€B(f,p)nF, ViG]0,e]. (5.17) 

In particular, there exists q\ £ {1, . . . , k} and i±, . . . , i qi £ I such that 

/ is strictly incoming to Hi q \fq = 1, . . . , q± (5.18) 

and 

f is parallel to H lq Vg > gi + 1. (5.19) 

Let F q = B(f , p)C\F with p > for g = 1, . . . , g x . Then p(F q ) > and it follows f rom fl5.18| ) 
that for p small enough, any / £ F q is strictly incoming to Hi q . Moreover, thanks to (5.17), any 
/ £ Fq is incoming to Hi for i £ J \ . . . , i qi }. Then we can use the induction hypothesis with 
x' = xo + e/o close to xo such that c(x' ) = k — qi < k to build F qi+ i, . . . ,Fk- The statement 
concerning the diameter of the Fj is a trivial consequence of the construction. □ 
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Corollary 5.4. Assume that (E,fj,) is weakly incoming to £1 and let xq £ Q. Then there exists 
r > 0, e > 0, F C E with fj,(F) > and 6 £ {±1} such that 

V/ £ F, VxG B(x , r ) n n, Vt£ [0, e], x + tOf £Q (5.20) 



Using these results and working as in Section 2 easily proves the following. 
Proposition 5.5. There exists N £ N and ci, C2 > suc/i that 

M h ( x , dy) = fih(x, dy) + cih~ d l\ x _ y \ <C2h dy 
where for all x £ Vt, fj,h(x, dy) is a positive Borel measure. 



(5.21) 



Proof. The starting point of the proof is to observe that for any k £ N and any non-negative 
function /, 

K k h f(x)> [ ... ! K h , ei ...K h>e J(x)d f i(e k )...d f i(e 1 ) (5.22) 



for any F\, . . . , F^ C E. Then the proof is the same as the proof of Proposition 2.2 In fact, (2.13) 



remains valid thanks to Lemma |5.1| Then we can mimick the end of the proof, using the fact that 
in Proposition 5.3 the set Fj can be chosen with arbitrary small diameter. Details are left to the 
reader. □ 



Proposition 5.5 implies a lemma analogous to Lemma |2.3| for the operator considered in 
this section. In particular, any function u £ L 2 (J7) satisfying 



\ 2 L2 + {{l-M h )u,u) L2 < 1 



(5.23) 



admits a decomposition u = ul + uh with < 1 and ||n//|| L 2 = 0(h). Using Proposition 5.3 



and the generalization of Lemma 2.3 gives Parts i and ii of Theorem 1 5. 2 



Part Hi is implied by the following lemma (where Bh still denotes the Dirichet form associated 
to M h ). 

Lemma 5.6. Let 9 £ C°°(r2) be fixed and let (tphi r h) £ H l (Q) x L 2 (f2) be such that \\rh\\L 2 (n) = 
O(h) and <~ph converges weakly in H 1 ^) to some (p. Then 



lim hr 2 B h (r h ,0) = 

/i->0+ 



and 



lim h 2 B h (f h ,8) = - [ [ V E f{e,x)V E 0{e,x) dxd^(e). 
h->o+ 6 J E J n 



(5.24) 
(5.25) 



Proof. The proof is the same as that of Lemma 2.4 



□ 



Total variation estimates for rates of convergence now follow as in Section 3 
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